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This  thesis  is  ray  meager  attempt  to  add  to  the  body  of 
knowledge  of  satellite  attitude  stability.  Due  to  the  immense 
cost  of  launching  satellites,  and  having  absolutely  no  place  on 
earth  to  test  stability  prior  to  launch,  only  satellites  which 
can  be  analyzed  for  stability  can  be  built.  I hope  that  perhaps 
this  study  may  be  of  use  to  those  who  design  satellites.  If  not, 
at  least,  I have  learned  a great  deal  about  satellite  stability 
and  its  problems. 
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This  study  considers  the  attitude  stability  of  a gyrostat 
satellite  containing  flexible  antennas,  in  which  any,  or  all, 
of  the  rotors  or  antennas  may  be  misaligned  with  the  main  body 
principal  axes.  The  problem  is  formulated  in  general  for  any 
number  of  rigid,  symmetrical,  spinning  rotors  which  are  fixed 
relative  to  the  body  and  for  any  number  of  antennas,  modeled 
as  rigid  rods  connected  by  torsional  springs.  The  stability 
analysis  is  based  on  the  Liapunov  direct  method,  using  the 
Hamiltonian  as  the  Liapunov  function.  Examples  are  presented 
for  a gravity-stabilized  satellite  containing  one  constant- 
speed  rotor  and  two  antennas  of  two  rods  each.  The  antenna 
misalignment  and  motion  is  restricted  to  the  xz  plane.  Results 
show  that  for  a rotor  misalignment,  the  largest  deviation  from 
the  position  defined  by  the  body  principal  axes  aligning  with 
the  orbital  axes  occurs  when  the  misalignment  is  toward  the 
orbit  tangent.  Misalignment  of  the  antennas  had  no  effect  on 
the  body  equilibrium  position  except  when  the  rotor  was  also 
misaligned.  In  that  case,  the  body  equilibrium  angles  were 
reduced . 
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ATTITUDE  STABILITY  OF  AN  ORBITING  SATELLITE 
CONTAINING  FLEXIBLE  ANTENNAS  AND  SPINNING  ROTORS 


I.  Introduction 


Background 

The  increasingly  complex  configurations  of  the  earth- 
orbiting satellites  have  generated  substantial  investigation 
of  satellite  attitude  stability.  Results  of  this  research 
provide  several  useful  analytical  methods  and  several  dif- 
ferent models  which  may  be  used  to  determine  stability. 

Pringle  (ref  1)  applied  Liapunov's  direct  method  to  determine 
the  stability  of  a satellite  with  connected  moving  parts. 
Meirovitch  and  Nelson  (ref  2)  used  an  infinitesmal  analysis 
on  a satellite  containing  elastic  parts  by  assuming  normal 
modes  of  elastic  displacement.  Nelson  and  Meirovitch  (ref  3) 
also  applied  Liapunov's  direct  method  to  a satellite  with 
elastically  connected  moving  parts  by  assuming  a lumped-para- 
meter representation  of  the  distributed  elastic  system. 

Likins  (ref  4)  made  a study  of  "dual-spin"  satellites  con- 
taining a spring-mass-damper  system  in  the  despun  portion  of 
the  satellite.  Crespo  da  Silva  (ref  5)  and  Rumiantsev  (ref  6) 
studied  satellites  containing  internal  rotors.  Rumiantsev's 
study  was  slightly  more  general;  however,  Crespo  da  Silva 
presented  results  which  were  both  more  extensive  and  in  a more 
useful  form.  More  recently,  Meirovitch  and  Calico  (ref  7) 
presented  an  elegant  method  to  determine  stability  of  a satellite 


containing  flexible  appendage::. . They  used  Liapunov's  direct 

method,  used  the  Hamiltonian  as  a Liapunov  functional,  and 
defined  integral  coordinates  to  evaluate  the  Hamiltonian. 

Calico  (ref  8,9)  extended  this  study  to  include  that  of  Crespo 
da  Silva  to  determine  stability  of  a satellite  containing  both 
flexible  appendages  and  internal  rotors.  These  references  are 
but  a few  of  the  numerous  studies  of  satellite  attitude  stabi- 
lity which  are  pertinant  to  this  study. 

Each  of  these  studies  has  one  common  characteristic! 
the  stability  was  investigated  about  a nominal  equilibrium 
position  defined  by  the  orbit  normal  and  the  body  angular 
momentum  vectors  being  aligned  with  a body  principal  moment- 
of-inertia  axis.  That  is  not  so  surprising,  for  generally, 
that  particular  position  is  the  one  of  primary  importance,  and 
if  it  is  not  stable,  the  satellite  is  considered  unstable.  In 
these  studies  it  was  generally  assumed  that  small  misalignment 
errors  in  the  rotors  or  antennas  do  not  affect  stability. 
However,  what  would  happen  to  stability  if  misalignment  were 
substantial,  e.g.,  an  extendable  antenna  failed  to  extend 
causing  a rotation  of  the  principal  moments  of  inertia? 

The  Problem 

The  problem,  then,  and  the  purpose  of  this  thesis  is  to 
determine  the  stability  of  a satellite  containing  both  flexible 
antennas  and  spinning  rotors,  any,  or  all,  of  which  could  be 
misaligned  with  the  body  principal  axes. 

The  General  Approach 

The  method  used  by  Calico  will  be  used  here,  i.e., 
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Liapunov's  direct  method  with  the  Hamiltonian  serving  as  the 

Liapunov  function.  The  potential  and  kinetic  energies  were 
determined  in  general  for  any  number  of  rotors  and  antennas, 
and  the  Hamiltonian  was  formulated  and  evaluated  for  sign 
definiteness  about  an  equilibrium  position.  For  the  specific 
example  presented,  the  equilibrium  positions  were  found  by 
setting  the  partial  derivatives  of  the  dynamic  potential  with 
respect  to  the  generalized  coordinates  equal  to  zero.  The  actual 
calculations  of  the  equilibrium  positions  and  the  evaluations 
of  the  Hamiltonian  were  done  numerically. 

Organization 

This  thesis  is  organized  so  that  long  mathematical  deriva- 
tions are  not  presented  in  the  text  but  are  included  as  an 
appendix.  Also  when  a large  set  of  equations  exists,  the 
entire  set  is  included  as  an  appendix. 


II . Mathcnnticnl  Formulation  of  the  Problem 


Model  and  Basic  Assumptions 

The  satellite  under  consideration  consists  of  a rigid  body 
containing  m rigid,  spinning,  symmetrical  rotors  and  has  attached 
to  it  n flexible  antennas.  The  spin  axes  of  the  rotors  are  fixed 
relative  to  the  main  body.  The  antennas  are  modeled  as  a number 
rigid  rods  connected  by  torsional  springs.  Damping  in  the 
satellite  is  not  modeled  explicitly?  however,  pervasive  damping 


Figure  1 . The  Satellite  Model 

will  be  assumed  for  those  configurations  which  are  otherwise 
stable,  thereby  insuring  assymptotic  stability. 

The  earth  is  modeled  as  the  central  body  of  mass  M 
producing  a spherically  symmetric  gravitational  force  field. 
The  total  mass  m of  the  satellite  is  small  compared  to  M. 

The  last  important  basic  assumption  is  that  the  motion 
of  the  center  of  mars:  of  the  satellite  is  unaffected  by  the 
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attitude  motion  or  the  elastic  motion.  This  assumption  is 

reasonable  if  the  dimensions  of  the  satellite  are  small  com- 
pared to  the  orbit  radius,  and  the  elastic  motion  is  either 
small  or  symmetrical  about  the  center  of  mass  C.  This  assump- 
tion allows  the  motion  of  C to  be  calculated  independently  from 
the  attitude  motion;  therefore,  the  orbital  motion  is  assumed 
to  be  known.  For  this  study,  the  satellite  is  in  a circular 


Fig.  2 depicts  the  coordinate  systems  used  in  this  thesis 
The  center  of  the  earth  is  assumed  as  the  origin  0 of  an 
inertial  reference  frame  XYZ.  The  center  of  mass  of  the  satel- 
lite in  its  undeformed  state  is  the  origin  of  a set  of  orbital 
axes  a^a^a^,  as  well  as  the  origin  of  the  body  fixed  axes  xyz . 
rij  is  ali'ned  along  the  position  vector  of  C in  XYZ;  a,,  is  the 
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orbit  tangent  vector  in  the  direction  of  the  velocity  vector; 

and  is  the  orbit  normal  vector.  When  the  satellite  is 

undeformed,  xyz  represents  the  principal  moment-of-inertia  axes. 

xiyizi  def^ne  the  aXes  for  "the  ith  antenna  and  are  fixes  with 

respect  to  xyz.  They  are  defined  such  that  when  the  antenna 

is  undeformed  it  lies  along  one  of  the  axes,  x.y.z.  define 

J J J 

the  nodal  axes  of  the  jth  rotor  such  that  z . is  both  the  sym- 

J 

metry  axis  and  the  spin  axis.  Because  x.y.z.  and  x.y.z.  are 

S 111  J J J 

fixed  in  xyz,  they  have  the  same  angular  velocity  as  xyz. 

Although  Fig.  2 shows  xyz  aligned  with  a^a^a^,  xyz  is  free  to 

rotate,  oi  denotes  the  origin  of  °j  denotes  the 

origin  of  x.y.z..  (Throughout  this  thesis  the  subscripts  "i" 

J J J 

denote  quantities  referred  to  the  antennas,  and  "j"  denote 
quantities  referred  to  the  rotors.  The  subscript  "o"  refers 
to  the  main  body.  "i"  numbers  from  1 to  n,  and  "j"  numbers 
from  n+1  to  n+m.) 

Fig.  3 depicts  the  position  vectors  and  their  components. 
Vector  quantities  are  indicated  by  an  underline,  and  their 
magnitudes  are  denoted  by  the  same  symbol  without  the  under- 
line, e.g.,  ^ is  the  position  vector  of  C in  XYZ,  and  Rc  is 
its  magnitude.  Time  derivatives  in  inertial  space  are  denoted 
by  a dot  above  the  quantity,  e.g.,  r is  the  time  derivative  of 
r in  inertial  space.  Time  derivatives  in  the  xyz  frame  are 
denoted  by  a dot  above  the  quantity  quantity  and  a prime  super- 
script, e.g.,  r'  is  the  time  derivative  of  r in  the  xyz  frame. 

The  System  Hamiltonian 

To  formulate  the  Hamiltonian  of  a dynamic  system,  the 

6 
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Figure  3.  Position  Vectors 


potential  and  kinetic  energies  must  be  determined.  The  poten- 
tial energy  will  be  determined  first.  Meirovitch  (ref  10*435) 
developes  the  gravitational  potential  for  a non-uniform  body 
in  a central  force  field.  To  a second  order  approximation 

= - GMrw  - GA^  [ (3f2-  l)(Iyy  * III  ~ If*  ) 

R HU 

+ (3ma-iyix*  +1*1  - Iyy)+  (3n*-l)^Iicx  +Iyy  “ In) 

+ I2(£mlxy  ♦ (n Iti  * vnrvlyi ) 

7 (l) 

v/here  1,  m,  and  n are  direction  cosines  between  xyz  and  , 

I , I , I ....  are  morpents  and  products  of  inertia  of  the 
xx'  yy  u K 
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ol'  the  satellite,  and  (J  is  the  universal  gravitational  constant. 


v 


The  total  gravitational  potential  of  the  satellite  is  the  sum  of 
the  potential  of  each  part.  Summing  over  the  body  and  rewriting 
eqn  (l)  in  matrix  form 

Vq  3 - GMvrt  - 36 H\  £.  Mt&ft 
3 fi*  Zjf? 


-fflfafij-l  ulUiuajfA.} 

where  Cii]  is  the  direction  cosine  matrix  between  311(1 

xyz , [j J is  the  inertia  matrix  of  the  ith  antenna  in  its 
deformed  state,  and  M is  the  direction  cosine  matrix 
between  and  xyz . The  elastic  potential  is  somewhat  easier 


to  write » 


Ve 


•itt  Kali-* 

£ 4«/  *•! 


(3) 


where  is  the  spring  constant  of  the  kth  spring  of  the  ith 

antenna,  and  |£ik  is  the  angle  through  which  the  spring  is 
displaced  from  its  undeformed  position.  Thus  the  total 
potential  is 

V * V3  + Ve  <»> 

The  kinetic  energy  can  be  written  as 

t 1 ij[  <5) 

where  m is  the  inertial  velocity  of  am  elemental  mass  dm, 
and  m is  the  mass  of  the  body.  When  m is  written  in  its 


. 
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component  terms , arid  the  internal  is  divided  into  i .•  Severn  i 

domains  DQ,  , and  the  kinetic  energy  becomes 

t-  * i MT/  UrUr«][A]M 

4*0 

♦ if  lartTillfl]-  4 £ I A M'lk'Uoi 


j*rvn 


<*7  -w 


(6) 


4*1 

+ WT  LtifCJjHAj] 

j*M» 

where  M is  the  velocity  of  Cj  /JL’l  is  the  mass  density  of 
the  ith  antenna » fo)  is  the  angular  velocity  of  the  jth 
rotor  relative  to  xyz;  [Uoi.  tPi)3  is  the  skew-symmetric  matrix 
which  yields  the  vector  cross  product  (r^+^Jxg?;  and  [j^-]  is 
the  inertia  matrix  of  the  jth  rotor  relative  to  x^y^^  Eqn 
(6)  is  derived  in  Appendix  A. 


Since  the  orbit  is  circular,  R_  is  constant,  and 

C 


sm  --  oo.1 

Rc 


(7) 


where  W.  is  the  orbital  angular  velocity.  The  angular  velocity 
of  the  xyz  frame  with  respect  to  ir  rtia  is 


[co]  “ fCoi]  + [u),] 


(8) 


where  is  the  direction  cosine  matrix  between  a^  and  xyz, 

and  M is  the  angular  velocity  between  xyz  and  a^a^a^* 


mwm 


Because  T and  V are  expi^ssed  in  a rotating  reference  frame 

and  the  orbital  velocity  is  constant,  the  system  is  considered 
non-natural,  and  it  is  explicitly  independent  of  time.  For  this 
system  the  Hamiltonian  assumes  the  form 


H-»  T, -T.  tv 


(9) 


where  T2  represents  the  terms  in  T which  are  quadratic  in  the 
generalized  veloxities,  and  TQ  represents  the  terms  in  T which 
are  independent  of  the  generalized  velocities.  Substituting 
eqn  (8)  into  eqn  (6)  and  grouping  terms 


To 


1*0 


(10) 


and 


Tz-  i 

ISO 

+ / Mi  lelffa'Jj  0i 

jsn+i  i*i  Joi 

+ Jl  [ Mi 


4T1  '0i 


(11) 


J*IVH 

Substituting  eqns  (4),  (10),  and  (11)  into  eqn  (9)  yields 
the  Hamiltonian  for  any  number  of  rotors  or  antennas. 

H-  H«.r£  UHtOUIM 

4*0 

+ i l {.Ail  * i t.  f m; 


(12) 


- ‘m 
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**  + [oj.n.  f A>  L(r*r+fr)]Ui]T{eJjo; 

C’l  y0i 

+ [tfiY  T*  L^]TLJi]{i^J 

- iwffiWf  LfcftaOWHM 

c-o 

-m  - trCUiTUMl) 

£«.  t*0 

- i Wo*  [Ha.,V(  f.  CfiaTt3ULCO)  {fa,] 

*i  £ £ Ku. n*u 
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Equilibrium  Position 

For  the  satellite  to  be  stable  about  an  equilibrium 
position,  H,  evaluated  at  that  position,  must  be  positive 
definite.  In  this  study,  because  the  rotors,  antennas,  or 
both  are  assumed  to  be  misaligned,  the  nominal  equilibrium 
position,  discussed  in  Chapter  I,  generally  does  not  exist. 
Therefore,  an  equilibrium  position,  hopefully,  one  near  the 
nominal  position,  must  be  found  and  H evaluated  for  sign 
definiteness.  The  Hamiltonian  may  be  expressed  as 

H = T,  + U 

where,  U,  the  dynamic  potential,  is  defined  as 

U ■ V-T,  <14 

Equilibrium  positions  are  defined  by  the  set  of  equations 

41/  - 0 d*i,n)  (15 

where  Xi  are  the  generalized  coordinates  and  n is  the 
number  of  coordinates.  In  general,  this  set  of  equations 
is  non-linear,  and  a numerical  solution  is  required.  Once 
the  equilibrium  position  is  found,  H may  be  evaluated  for 
positive  definiteness  to  determine  stability. 

Evaluation  of  the  Hamiltonian 

For  H to  be  positive  definite,  only  U need  be  evaluated 
because  by  definition  T^  is  positive  definite;  therefore, 


1.  +*  ^ 


is  dele rained  by 


H is  greater  that  U.  Ji.  n definiteness  ci' 

forming  the  matrix  [u^] , where 

U i 

and  applying  Sylvester's  criterion,  i.e.,  for  positive  definite- 
ness the  determinant  of  each  of  the  leading  principal  minors 
of  must  be  greater  than  zero.  If  U is  positive  definite, 

then  the  equilibrium  position  is  stable,  and,  moreover,  it  is 
assymptotically  stable  when  pervasive  damping  exists. 


£li  I . . 

'dh^Lj  (equilibrium 


IV.  Stability  of  Specific  Examrl.es 
Constant  Speed  Rotor 

In  each  of  the  specific  examples,  the  rotor  is  assumed  to 

be  driven  at  a constant  speed.  Since  X.y.z.  are  nodal  axes, 

J J J 

= (o  0 s -1 T where  s-  is  the  constant  spin  rate  of  the  rotor. 

J J 

When  this  is  the  case,  the  last  term  in  T2  reduces  to  a Tj  term 

and  a Tq  term.  Since  neither  Tj  nor  T2  appear  in  U,  only  the 

T term  will  be  considered.  The  additional  T term  becomes 
o o 


<Aff<u3TI  Ui3rt3j{ni3 


j 

Therefore,  Tq  for  a constant  speed  rotor  becomes 

Toe  - To  two  LTjI 

j=nti 

For  all  calculations  TQC  will  be  used  for  Tq  in  eqn  (14). 


(17) 


Examples 

Three  specific  examples  are  considered.  The  first  is  a 

satellite  containing  one  constant  speed  rotor  whose  spin  axis  is 

nominally  aligned  with  the  z-axis,  but  misalignment  is  allowed. 

% 

The  second  example  contains  one  constant  speed  rotor  aligned  with 
the  z-axis,  but  it  has  tv/o  antennas  restricted  to  symmetrical 
motion  but  which  may  be  misaligned.  The  last  example  is  a com- 
bination of  the  first  two.  Each  antenna  is  modeled  by  two  rods 
of  length  L,  mass  m/2,  and  is  connected  by  two  springs  of  stiff- 
ness K.  The  antennas  are  attached  to  the  main  body  at  a distance 
♦h,,  in  the  x-direction  and  h in  the  z-direction.  The  antennas 
are  further  restricted  to  motion  only  in  the  xz-plane.  The 
model  is  shown  in  Fig.  B-l  in  Appendix  B. 
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For  convenience  of  calculation,  the  matrices  will  bo 

redefined  as  and  where  r denotes  the  undeformed 

state,  and  e denotes  the  change  due  to  elastic  displacements, 
r^r  and  [Jje  are  developed  in  Appendix  B.  ty  r will  then 
be  summed  with  cy  to  form  [jJr,  the  inertia  matrix  of  the 
entire  undeformed  body. 


[J]T  1 [«  +£  ILTltirltil 


By  appropriately  defining  the  principal  moments-of-inertia 


/»•  O 

fa]  = -.0* 

^ c. 


- 

fUrMrUl-mL*  ^ Xlr° 

L O XjrJ 


where 


I.r  * 2( 

Ur  - 2(^,4^  1 
hr  s I \r  + lar 

The  inertia  matrix  of  the  undeformed  body  becomes 

^•  + mLaI*c  O /l  O 

[j]r  = 0o+mllIar  - 0 

, ^ Co  -V  Ml.1 1 jr.  .O  C 


r 


Also 


[j]«  = ml?  L* hP 

O Iae 


where 


L« 3 -3  sty  + $ s'ty  +^sl0i  - r S4>x  + s0,s0z 

lie  5 3^? C0,+ + ^ C&+C0,C&-  § (24) 

Xae  * Xi«  + lae. 

For  the  single  rotor»  define  [j~l  and  iXijj. 

r 1 °1 

[jj]  = I (25) 

[o  jJ 

fa»r'{ooj3  <26) 

With  the  inertia  matrices  defined,  and  the  rotor  speed  defined, 
the  direction  cosine  matrices  must  be  expressed  explicitly 
to  enable  U to  be  calculated. 

For  a set  of  Euler  angles  defined  by  a 02  rotation  about 
y,  a -0^  rotation  about  x,  and  a 0^  rotation  about  z,  and 

(M  become 

C0,c©j  “ i0,5©i5ej  [ 

f Ca.l  - - (cexse3  + s6lsexcej)  f , . 


(U 


ce»5©*. 

'(se*ce3  4 se,c©4s©j 

SOzSOj  “ 50,00x003 

ce.ce>. 


Also 


i a]  s m-  fxK] 

iu-C-8,] 


Colt 

Li  j]  = -CA.SUL 


SoCt 

O 

CcL\  Cdj, 

-<S*» 

S^CoCj, 

<Jo(i 

L-<Sx,^o(j.  MjCoCi,  <J oC%  J 

where  a<2  is  a rotation  about  and  is  a rotation  about 

Xj. 

Substituting  eqns  (2) , (3) , (7) , (10) , and  (1?)  <nto  eqn  (14), 
dividing  by  Cul£ , and  performing  the  indicated  matrix  operations 
yields  the  scalar  equation  for  UQ. 

Ucs  +lL(ceaceJ-s©,^5©i)lCA*^i5ilifi) 

+(ce>se j*  se,sejcej)l(Qr&gja9)  + (ce, se*)z(c + m lz lie )] 

+ i [-  (S6*ce*  + se, C61  j0i)V>U  wLlI.^ 

c 

- - 5€\c©ace,)a(  a-»ml3Iac^  (29) 

c 

- ( ce, c©Qa( C4  ml5i 3e)l - p[(s©,cej  4. 

o 

+ ^6>aS©j  - 50»C©i c©i)*S?(  1 Cc<i  + C0,C0aCo<i} 

4 wlMi;  +4i) 
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program  was  designed  to  evaluate  the  equilibrium  equations  and 
to  test  the  sign  definiteness  of  The  equilibrium  equa- 

tions  and  LL  ^ are  listed  in  Appendix  C. 

In  addition  to  0,  A AOq*  . and  <<2  • following  para- 
meters were  defined. 

Kj  = (C-B)/A 

K2=(C-A)/B 

HX=hx/L 

HZ=hz/L 

RA=mL2I3r/Ao 

2 2 

Stability  was  tested  for  various  values  of  , oC^ , 0,  A A>0» 

HX,  HZ,  RA,  , and  and  K2  were  varied  throughout  the 

region  indicated  in  Fig.  4 by  the  diagonal  lines.  Crespo  da 
Silva  and  Calico  found  the  regions  of  stability  for  the  nominal 
equilibrium  position  to  lie  within  this  area. 


V.  Result;;  and  Cone  1 u::  ions 


Results 

Obviously  there  are  numerous  solutions  to  the  equilibrium 
equations;  however,  only  the  stable  equilibrium  positions  near 
the  nominal  position  are  considered. 

Example  I.  For  the  first  example  RA=0;  therefore,  the 
only  parameters  which  were  varied  were  , K2 , P,  , and  <x^. 

When  2~°  • i»e.,  when  there  is  no  misalignment  of  the  rotor, 
the  results  duplicate  those  of  Crespo  da  Silva  for  the  deriva- 
tion is  identical.  However,  when  0^  does  not  equal  zero,  the 
nominal  equilibrium  position  no  longer  exists.  It  is  interesting 
to  note  that  the  stable  equilibrium  is  still  defined  by  the 
angular  momentum  vector  and  the  orbit  normal  vector  aligning! 
however,  they  are  no  longer  aligned  with  a body  principal  axis. 
This  is  easily  demonstrated  for  the  simple  case  where  tf2=-90°. 

For  this  case  the  equilibrium  equation  reduces  to 


S02C&X  + 0(c vifiiCol,)  * o 

c 

where  ©^=0^=0.  The  momentum  is  defined  by 

[l]  - LJ]r  [h*]  tOo  * 


(30) 


(31  ) 


which  reduces  to 


r, 


(4  = 


(-  60e i&xCG*  *tdoCj6aC©.a  t 3s  Soi ; C&j  + Js  CoLi  JSa') 


U)e  /) s,x6x  ' U) o x + 3i  C*t  COx  + 3s  S*,  sex) 
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(32) 


f 


The  first  terra  in  eqr,  is  exactly  t p.  v>  , , 


l *_A. 


u . ar 


momentum  of  the  equilibrium  position  is  aligned  with  the  orbit 
normal . 


Fig.  5»  6,  7,  and  8 show  stability  regions  and  equilibrium 
angle  curves  in  the  plane  for  various  values  of  , otg  and 

For  given  values  of  and  3,  equilibrium  angles  furthest 
from  the  nominal  were  found  to  occur  when  o<2=0,  or  "the  misalign- 
ment of  the  rotor  was  toward  the  a2  axis.  When  misalignment 
was  toward  the  a1  axis,  the  equilibrium  position  is  very  near 
the  nominal.  Fig.  5.  6,  and  8 also  show  that  ®2e=_^l  811(1  ®le=0 
when  K^=0  or  K2=l . Fig  7 shows  that  0^=-*^  and  ®2e=®  when 
K2=0.  In  both  cases  9^=0.  In  fig.  5«  6.  7,  and  8 the  region 
above  the  solid  line  K^Kg  in  the  first  quadrant  and  above  the 
solid  curve  in  the  second  quadrant  is  the  stability  region  for 
the  nominal  equilibrium  position.  In.  fig.  5 and  6,  the  region 
below  the  solid  curve  in  the  second  quadrant  is  still  stable 
about  an  equilibrium  position,  but  the  equilibrium  angles  be- 
come too  large  to  be  of  interest.  A comparison  of  fig.  5 and 
6 shows  that  there  is  essentially  no  difference  between  the  two 
for  areas  I and  II;  however  for  0^=10°,  area  III  is  very  slightly 
larger  than  that  for  0^=1°.  The  boundary  curve  of  area  III  in 
fir*.  6 crosses  the  Kj^  axis  and  Kj=-1  line  at  nearly  the  same 
place  as  in  fig.  5s  however  the  curve  is  slightly  flatter. 

This  difference  is  a little  more  obvious  in  fig.  9 where  0. 

1 e 

is  plotted  vs.  -0^  for  a give  value  of  K.  , 3,  and  , and  four 

different  values  of  Kg,  e.g.,  for  K2=.2,  -0^=1°,  and 

-<X,  r 10°,  9 -29°.  For  small  values  of  -ol  , 6,  is  nearly  a 

1 1 e lie 
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ror  § 1,  the  re;  ion  of 


} 


! 


linear  fur.ction  of  Kr  any 

stability  is  considerably  larger  than  for  $-.2.  With  much  more 
angular  momentum  in  the  rotor,  the  equilibrium  angles  are  very 
much  closer  to  the  misalignment  angle,  -0^ . 

Presentation  of  data  for  misalignment  when  o<2  is  different 
from  zero  or  -90°  is  difficult  because,  in  this  case,  all 
three  Euler  angles  are  different  from  zero.  However,  for  all 
stable  equilibrium  positions  which  were  found,  0^e  was  aJhmys 
less  than  .01°.  Analysis  of  the  numerical  data  suggest.  $hat 

©le  is  approximately  equal  to  -a^cosWg  when  ^=0  or  K2^4*»  and 
02e  is  approximately  equal  to  o^sind^  when  1^=0.  To  confirm 
this  observation,  it  is  useful  to  consider  two  other  relation- 
ships between  and  K2,  and  the  equilibrium  equations. 

= Ka(Kv-J) 

C K,Ki-l 

c-a  - k>(K2~0 

C K»Ka-| 


(33) 

(34) 


(&£)lCGtS26tSZ6i  - <M5let-/)  ^©j] 

- (c-$)  -i)  - 

c 

- CGi  C©*  C©i  Sdi  CoU  - 50,  ce*  c = o 


(35) 


2(gt&)  Cs©,c20a520j  + wezCc'ej-s’e.s2©,)] 
c 

+2(Qs&) C1©, 52©z  -^CCCQaC©,  - J©.S©tS©j)  MyScti 

C 

(C©aS6j  4 3©i  S©i  - CG,  56*  C*,*]  ■*  O 
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( y ) 


■ 


(37) 


Let  e-je-°  and  (C- B)/C- 0 , then  eqn  (35)  becomes 

C0,C0t  * S©,C©*Co(,  = O 

which,  for  small  and  9le>  reduces  to 

Qie  ~ -o(, Cola.  (38) 

Also,  letting  e.je=0,  (C-A)/C=0,  and  ©le  be  small,  eqn  (36) 
becomes 

C0a  S*,Scit-  CQiSOx  Crf,  =0  (39) 

which  reduces  to 

02ft  * (40) 

Therefore,  when  p(2  is  different  from  zero  or  -90°  and  $ = .2, 

©le  is  determined  from  fig.  5 or  6 , and  e2g  is  determined 
from  fig.  7s  is  assumed  approximately  equal  to  zero.  To 

determine  ©le  from  fig.  5 or  6 , replace  by  e^coscfc,  in  the 
legend  for  I,  II,  and  III.  To  determine  ©2e  from  fig.  7, 
replace  by  -at^sinoig.  These  approximations  are  within 

five  percent  of  the  numerical  data. 

Example  II.  The  second  case  considered  a satellite  with 
no  rotor  misalignment  but  a misalignment  HZ  in  the  z-direction. 
The  specific  purpose  of  this  example  was  to  determine  if  a 
misalignment  out  of  the  xy-plane  could  cause  instability. 

For  this  example  HZ  was  varied  from  zero  to  one  for  HX=.2, 

0 .2,  KA^.25,  and  X ? /<0?Q= . 11 33-  Cases  for  HX^.5,  RA  .5,  and 
/o;^- • 2?f’7  were  also  studied.  Table  1 shows  that  for  HZ- 0 
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4 the  satellite  was  stable-  about  the  nominal  position  where  all 

angles  were  equal  to  zero.  As  HZ  was  increased  to  one,  and 

gig  increased  significantly;  however,  6le»  ®2e  ’ 5111(1  03e  no * 

change  from  zero,  i.e.,  the  body  orientation  remained  the  same, 
but  the  antennas  are  no  longer  stable  in  an  equilibrium  posi- 
tion of  zero  energy  in  the  springs.  For  HZ/O , 0^  and  0^  are 
found  in  Table  I,  and  the  stability  region  for  the  body  is 
found  in  fig.  5 (the  nominal  stability  region). 


Table  I. 

Ranges  of  0^  and  0, 

) 

*1 

*2 

HX 

HZ 

RA 

1 

2 

0 

0 

.2 

0 

.25 

.1133 

0 

0 

0 

0 

.2 

.1 

.25 

.1133 

4.74 

1.09 

0 

0 

.2 

• 5 

.25 

.1133 

23.74 

5.99 

0 

0 

.2 

1 .0 

.25 

.1133 

46.88 

16.13 

0 

0 

.2 

.5 

.5 

.1133 

23.75 

6.00 

5 

0 

.2 

0 

.25 

.1133 

0 

0 

5 

0 

.2 

.1 

.25 

.1133 

4.74—4.53 

1 .08—1 .19 

5 

0 

.2 

• 5 

.25 

.1133 

23.74-23.08 

6.00—6.42 

5 

0 

.2 

1 .0 

.25 

.1133 

46.87-45.13 

16. 13-16. 46 

5 

0 

.2 

.5 

• 5 

.1133 

23.74-22.90 

6.00—6.40 

5 

0 

.2 

.5 

.5 

.2267 

21.90-21 .00 

8.13—8.47 

5 

5 

.2 

0 

• 25 

• 1133 

0 

0 

5 

5 

.2 

.1 

.25 

.1133 

4.75—4.53 

1 .09--1 .19 

5 

5 

.2 

• 5 

.25 

.1133 

23.74-22.73 

5 . 99--6 .42 

5 

5 

.2 

1.0 

.25 

.1133 

46.87-45.17 

16. 13-16. 46 

5 

5 

.2 

.5 

.5 

• 1133 

23.47-22.70 

6.00—6.51 

10 

0 

.2 

.5 

.5 

.2267 

21.90-20.70 

8.18—8.54 

t* -2 

for 

all 

values  in  this 

table 

Example  III.  As  shown  in  Table  I,  when  HZ-0,  0^-0^=O 
is  still  a stable  equilibrium  position.  The  numerical  data 
also  showed  that  the  body  equilibrium  angles  did  not  change  from 
those  without  the  antennas.  Therefore,  fig.  5i  6,  7.  and  8 
ern  ve  used  to  determine  stability  regions  and  equilibrium 
• 1 <•»  curves  for  th'  satellite  with  or  without  antennas 
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provided  HZ-0.  When  HZ/O,  the  body  equilibrium  ancles  were 

reduced.  As  HZ  was  increased  the  stable  equilibrium  angle 
curves  shifted  to  the  left  as  indicated  in  fig.  10.  Also  as 
RA  was  increased  to  .5,  the  equilibrium  angle  curves  again 
shifted  slightly  to  the  left  as  indicated  in  fig.  11.  Doubling 
the  value  of  A /too  did  not  have  significant  effect  on  the 
body  equilibrium  angles,  but  was  reduced  by  about  8%,  and 
02  was  increased  by  about  30%.  and  $2  were  very  weak  func- 
tions of  and  K2,  and  Table  I shows  the  range  through  which 
^ and  02  varied  as  ^ and  K2  were  varied. 

Conclusions 

Fig.  5,  6,  and  8 clearly  indicate  that  if  the  satellite  is 
designed  for  stability  about  the  major  principal  axis,  rotor 
misalignment  has  much  less  effect  on  the  equilibrium  angle  than 
if  designed  for  the  region  where  is  less  than  zero.  In  this 
case,  the  equilibrium  angle  is  always  less  than  the  misalignment 
angle.  If  stability  is  designed  for  the  region  between  the 
nominal  equilibrium  curve  and  region  III,  very  minor  misalign- 
ment of  the  rotor  can  cause  significant  changes  in  the  equil- 
ibrium angles.  For  misalignment  in  the  antennas,  the  conclu- 
sions are  not  so  clearly  defined  because  the  specific  model  was 
very  restricted.  For  this  model,  the  misalignment  actually 
reduced  the  body  equilibrium  angles  if  there  was  a rotor 
misalignment,  and  had  no  effect  if  there  was  no  rotor  misalign- 
ment. It  was  tacitly  assumed  that  the  antenna  equilibrium 
an  les  wore  not  important  to  the  actual  use  of  the  antenna. 

?7 


■ '*vW 


fir 


Figure  10.  Body  Equilibrium  Angle  Comparison  for 

Changes  in  HZ  •> 

* i=5°.  <X2=0.  $ = .2,  HX= . 5 , RA=.25,  A/k‘=.U33 
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The  theory  presented  is  valid  for  more  General,  less  restricted 
antenna  models.  For  instance,  misalignments  other  than  HZ, 
asymmetric  motion,  or  motion  in  the  xy-plane  could  be  allowed. 

To  do  this,  one  would  only  have  to  redefine  the  inertia  matrices 
of  the  antennas . 


29 


i M i o rn  nhv 


1.  Pringle,  R.  Jr.,  "On  the  Stability  of  a Body  with  Connected 
Moving  Parts,"  AIAA  Journal . Vol . 4,  No.  8,  1966,  1 395-1 404. 

2.  Meirovitch,  L.  and  Nelson,  H.D.,  "On  the  High-Spin  Motion 
of  a Satellite  Containing  Elastic  Parts,"  Journal  of  Space- 
craft and  Rockets.  Vol.  3,  No.  11,  1966,  1597-1602. 

3.  Nelson,  H.D.  and  Meirovitch,  L. , "Stability  of  a Nonsym- 
metrical  Satellite  with  Elastically  Connected  Moving  Parts," 
The  Journal  of  the  Astronautical  Sciences,  Vol.  13,  No.  6, 
226-23^- 

4.  Likins,  P.W.,  "Attitude  Stability  of  Dual-Spin  Systems," 
Hughes  Aircraft  Company.  Space  Systems  Divisions  Report. 

SSD  60377R,  September  1966. 

5.  Crespo  da  Silva,  M.R.M.,  "Attitude  Stability  of  a Gravity- 
Stabilized  Gyrostat  Satellite,"  Celestial  Mechanics.  Vol. 

2,  1970,  147-165. 

6.  Rumiantsev,  V.V.,  "On  the  Stability  of  Stationary  Motions 
of  the  Gyrostat  Satellite,"  18th  Congress  of  the  Inter- 
national Astronautical  Federation.  Belgrade.  1967. 

7.  Meirovitch,  L.  and  Calico,  R.A.,  "The  Stability  of  Motion 
of  Satellites  with  Flexible  Appendages,"  NASA  Contractor 
Report . NASA  CR-1978,  February  1972. 

8.  Calico,  R.A.,  "Stability  of  a Gyrostat  Satellite  with 
Flexible  Appendages,"  Journal  of  Spacecraft  and  Rockets. 

Vol.  13,  No.  8,  August  1976,  505-508. 

9.  Calico,  R.A.,  "Stability  of  a Flexible  Satellite  Containing 
Internal  Rotors,"  AIAA  Paper  76-786 . AIAA/AA3  Astrodynamics 
Conference,  San  Diego,  California,  August  1976. 

10.  Meirovitch,  L. , Methods  of  Analytical  Dynamics.  McGraw-Hill 
Book  Co,  N.Y. , 1970. 


30 


99QJIPP 


Appendix  A 


Kinetic  Energy  Derivation 

The  kinetic  energy  of  a body  may  be  written  as 

T - if  {RfCft} <5> 

*fV\ 

If  then  eqn  (5)  becomes 

T -if" * i[jrV{rUm  ♦ [ftffjrUm  <41  > 

By  definition  of  the  center  of  mass  L frlclm  3 0 , and  the  last 
term  in  eqn  (41)  is  identically  zero.  There  are  several 
distinct  domains  in  the  body,  and  the  second  term  in  eqn  (41 ) 
is  more  easily  calculated  by  integrating  over  each  domain  and 
summing  over  the  entire  body. 

Considering  the  rigid  portion  first,  r = <*>  x rQ  because 
r*  = 0.  Therefore 


i l [r]TMcU#  = j {r0]r[co}  <U  fwj 

■ i MfWM 


(42) 


So 


and 


For  the  ith  antenna,  r = u>  x r.  + r.' ; however  r.  = r . + o.  . 

l *“1  —I  — 01  r"l 

r * &*<ui  * i[  (43) 

ij*.  fr3T{r3cU;  si  (wj rJm.%  [fc]T[C«*;  +^)]  [(r*i  *?;.)][ 


(44) 
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The  l'irst  term  on  the  right  of  the  equal  sign  in  eqn  (44) 

reduces  to  ^ £w3  [£‘il 

Similarly  for  the  jth  rotor,  r = w x r'. , and 

J 

x ly  So 

i L-  fr]Tif3<Uj  i toru]iaj*MWuat«M  <45> 

The  part  of  the  kinetic  energy  of  the  rotor  which  does  not 
arise  from  its  spin  has  been  included  in  [JQ3 . 

Combining  eqns  (5)»  (42),  (44),  and  (45)  yields  the  total 
kinetic  energy  for  the  system  (eqn  (6)). 
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Apt i.  ri'lix  b 


Formulation  of  [j/]r  andfjjlg. 

Fig.  B-l  shows  the  model  used  for  the  antennas.  Since  they 
are  symmetrical  about  the  z-axis,  and  are  restricted  to  motion 
in  the  xz-plane,  all  products  of  inertia  vanish,  and  I is  the 
sum  of  Ixx  and  T-zz>  Since  the  antennas  are  symmetrical,  only  one 


Figure  B-l . Antenna  model 

is  shown.  From  the  definition  of  moments-of- inertia 
Ixx=2jIy<z'>dD;  therefore, 

fL  r2L 

I**  5 2J  +2J  2:s4>L  - Ls<t,)  dt 

■'o  ■'J. 


~ hft  La(H2*  t sx<t>,  t -HZi&i  +s<$'S0i) 


(46) 
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(4?) 


2JlmO u + + Lc<t,)d* 

i 2 H X c ■♦•HXcCJi  + C$,C0t 


[ J.V  consists  of  those  values  of  I , I , and  I evaluated 
for  <&^  = = 0.  fjjg  consists  of  the  complete  matrix  £l^ 

less  which  yields  eqns  (21)  and  (24).  It  should  be  noted 

this  inertia  matrix  for  the  antennas  is  valid  only  for  symmetri- 
cal motion. 


Equilibrium  Equations  and  Dynamic  Potential 


Equilibrium  Equations 

The  following  equations  were  used  to  determine  the 
equilibrium  angles. 


-(c'-3')s|ei  ( vs^x-i)  - p[e©,ce*seti 


(48) 

i 


- ce,C©iCea  ,S<*»CoU  - 50lC6a.Cot,3  = Q 


&L  - Z [ se,  do*  s ze»  4 s&i  (c les  -s1q,61q4 )] 

+ cle,  s^ea  - pUceicej-se^easej).^#,.  (49) 
c 

4(cej5©a  4 ss.ieace^s^c**  - ce,  5e>co<i3  - 0 


dJl.  - (a^;)  [2se,s2fl,.C20a  4 szejOlcle>'i) 
Mi  c 2 

-5f©.  5aejC^5*©»-/)]  - ^(-JeaS©^  + 

2 


4-  (S©aC©4  + S©,c©as©» ).*<*,  C^J  = O 


(50) 


r 


= 1 (C0xCGj-i©^©i6ej)l>  as'  (coxsej 
2C  a<p, 

f j©,s©%c©,)1  + (cetset)*y  j^lCsetcej 
<3«i 

4-  5e.ce**©,)1  ♦ £b'  (seasaj-s^coxc^i)2,  (51) 

64,  64, 
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Dynamic  Potential 

The  following  equations  are  the  dynamic  potential  used  to 
formulate  and  determine  its  sign  definiteness. 
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